home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_asum.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  2.8 KB  |  144 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_asum (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int N = 1;
  13.    float X[] = { 0.239f };
  14.    int incX = -1;
  15.    float expected = 0.0f;
  16.    float f;
  17.    f = cblas_sasum(N, X, incX);
  18.    gsl_test_rel(f, expected, flteps, "sasum(case 40)");
  19.   };
  20.  
  21.  
  22.   {
  23.    int N = 1;
  24.    double X[] = { -0.413 };
  25.    int incX = -1;
  26.    double expected = 0;
  27.    double f;
  28.    f = cblas_dasum(N, X, incX);
  29.    gsl_test_rel(f, expected, dbleps, "dasum(case 41)");
  30.   };
  31.  
  32.  
  33.   {
  34.    int N = 1;
  35.    float X[] = { 0.1f, 0.017f };
  36.    int incX = -1;
  37.    float expected = 0.0f;
  38.    float f;
  39.    f = cblas_scasum(N, X, incX);
  40.    gsl_test_rel(f, expected, flteps, "scasum(case 42)");
  41.   };
  42.  
  43.  
  44.   {
  45.    int N = 1;
  46.    double X[] = { -0.651, 0.079 };
  47.    int incX = -1;
  48.    double expected = 0;
  49.    double f;
  50.    f = cblas_dzasum(N, X, incX);
  51.    gsl_test_rel(f, expected, dbleps, "dzasum(case 43)");
  52.   };
  53.  
  54.  
  55.   {
  56.    int N = 2;
  57.    float X[] = { 0.899f, -0.72f };
  58.    int incX = 1;
  59.    float expected = 1.619f;
  60.    float f;
  61.    f = cblas_sasum(N, X, incX);
  62.    gsl_test_rel(f, expected, flteps, "sasum(case 44)");
  63.   };
  64.  
  65.  
  66.   {
  67.    int N = 2;
  68.    double X[] = { 0.271, -0.012 };
  69.    int incX = 1;
  70.    double expected = 0.283;
  71.    double f;
  72.    f = cblas_dasum(N, X, incX);
  73.    gsl_test_rel(f, expected, dbleps, "dasum(case 45)");
  74.   };
  75.  
  76.  
  77.   {
  78.    int N = 2;
  79.    float X[] = { -0.567f, -0.645f, 0.098f, 0.256f };
  80.    int incX = 1;
  81.    float expected = 1.566f;
  82.    float f;
  83.    f = cblas_scasum(N, X, incX);
  84.    gsl_test_rel(f, expected, flteps, "scasum(case 46)");
  85.   };
  86.  
  87.  
  88.   {
  89.    int N = 2;
  90.    double X[] = { -0.046, -0.671, -0.323, 0.785 };
  91.    int incX = 1;
  92.    double expected = 1.825;
  93.    double f;
  94.    f = cblas_dzasum(N, X, incX);
  95.    gsl_test_rel(f, expected, dbleps, "dzasum(case 47)");
  96.   };
  97.  
  98.  
  99.   {
  100.    int N = 2;
  101.    float X[] = { 0.169f, 0.833f };
  102.    int incX = -1;
  103.    float expected = 0.0f;
  104.    float f;
  105.    f = cblas_sasum(N, X, incX);
  106.    gsl_test_rel(f, expected, flteps, "sasum(case 48)");
  107.   };
  108.  
  109.  
  110.   {
  111.    int N = 2;
  112.    double X[] = { -0.586, -0.486 };
  113.    int incX = -1;
  114.    double expected = 0;
  115.    double f;
  116.    f = cblas_dasum(N, X, incX);
  117.    gsl_test_rel(f, expected, dbleps, "dasum(case 49)");
  118.   };
  119.  
  120.  
  121.   {
  122.    int N = 2;
  123.    float X[] = { -0.314f, -0.318f, -0.835f, -0.807f };
  124.    int incX = -1;
  125.    float expected = 0.0f;
  126.    float f;
  127.    f = cblas_scasum(N, X, incX);
  128.    gsl_test_rel(f, expected, flteps, "scasum(case 50)");
  129.   };
  130.  
  131.  
  132.   {
  133.    int N = 2;
  134.    double X[] = { -0.927, 0.152, -0.554, -0.844 };
  135.    int incX = -1;
  136.    double expected = 0;
  137.    double f;
  138.    f = cblas_dzasum(N, X, incX);
  139.    gsl_test_rel(f, expected, dbleps, "dzasum(case 51)");
  140.   };
  141.  
  142.  
  143. }
  144.